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Abstract 



We present a continuous time quantum search algorithm analogous to Grover's. 
In particular, the optimal search time for this algorithm is proportional to y/~N, 
where ./V is the database size. This search algorithm can be implemented using 
any Hamiltonian with a discrete energy spectrum through excitation of resonances 
between an initial and the searched state. This algorithm is robust and, as in the 
case of Grover's, it allows for an error 0(l/yN) in the determination of the searched 
state. A discrete time version of this continuous time search algorithm is built, and 
the connection between the search algorithms with discrete and continuous times is 
established. 
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1 Introduction 

In the last twenty years the attention of researchers of several different areas 
has been attracted towards quantum computation [1,2]. This field of knowl- 
edge presents new scientific challenges both from the theoretical and the ex- 
perimental points of view. Part of the theoretical challenge is to learn how to 
work with quantum properties to obtain new and more efficient algorithms. 

A quantum algorithm can be seen as a definite sequence of unitary transfor- 
mations acting over a quantum state, in some Hilbert space. Its size is pro- 
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portional to the number of elementary unitary transformations of which the 
algorithm is composed. Concepts like interference phenomena, quantum mea- 
surements, resonances, quantum parallelism, amplification techniques, etc., are 
employed by this new computation science. However, relatively few quantum 
algorithms were created; among them, Shor's and Grover's [3,4] algorithms 
are the best known. Shor's algorithm decomposes a number in its prime fac- 
tors more efficiently than any known classical algorithm. To achieve this it 
uses quantum parallelism, quantum Fourier transforms, and the properties of 
quantum measurements. Grover's search algorithm locates a marked item in 
an unsorted list of N elements in a number of steps proportional to y/N, in- 
stead of O(N) as in the classical case. It performs a unitary transformation 
of the quantum state which increases the likelihood that the marked state 
of interest will be measured at the output (amplification technique). It has 
been proven that there are neither quantum nor classical algorithms that can 
perform faster such an unstructured search [5]. 

In this work we present a continuous time quantum search algorithm, which is 
controlled by a time dependent Hamiltonian. In particular there are not uni- 
tary operators that are applied at certain time steps, as in Grover's algorithm. 
Unlike other authors [6,7] who studied continuous time search algorithms, the 
most relevant characteristic of our model is the use of quantum resonances, 
showing explicitly Grover's assertion that his algorithm is a resonance phe- 
nomenon [8]. We should remark that a suggestion to employ resonances in 
quantum computing was made by M. Zak in [9]. The paper is organized as 
follows. In the next section we develop the continuous time search model. In 
section 3 we present examples of use of this model for two typical Hamiltoni- 
ans. Then we show, in section 4, the connection of this continuous time with a 
discrete time search algorithm. In the last section, we make some concluding 
remarks. 



2 Resonances 

The evolution of the wavefunction \^(t)) satisfies the Schrodinger equation, 

where H = H + V(t) and we have taken Planck's constant h — 1. Here 
H is a known nondegenerate time-independent Hamiltonian with a discrete 
energy spectrum, and V(t) is a time dependent potential that will be defined 
below. We should point out that the extension to the degenerate case is quite 
straightforward, as in the case of Grover's algorithm for several equally marked 
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items. 



Let us consider the normalized eigenstates {\n)} and eigenvalues {e n } for the 
unperturbed Hamiltonian. These sets can be finite or infinite, depending on 
the nature of H . We now consider a subset N of {\n)} formed by N elements 
on which we shall apply the search algorithm. We take a known eigenstate \ j), 
with eigenvalue Ej, as the initial state of the system. This initial state does 
not belong to the search set N. 

Let us call \s) the unknown searched state whose energy e s is known. This 
knowledge is equivalent to 'mark' the searched state in Grover's algorithm. 
Our task is to find the eigenvector \s) which transition energy from a given 
initial state \j) is cu S j = Ej — e s . We propose the following potential V, which, 
as it may be easily verified, produces a resonance between the initial and the 
searched states, 



V(t) = \p) (j\ exp (iu sj t) + \ j) (p\ exp (-iu sj t) , (2) 
where \p) = ~^'^2\ n ) is an unitary vector, that can be interpreted as the 

rteN 

average of the set of vectors N. This definition assures that the interaction 
potential V is hermitian, that the transition probabilities from state \ j) to any 
state of the set N are all equal, W n j = \(n\V(t)\j)\ 2 = l/N, and finally that 
the sum of the transition probabilities verifies ^ W n j = 1. 

Let us express as an expansion in the eigenstates {\n)} of H , \^f(t)) = 

J2m a m{t) exp (— ie m t) \m), where the expansion coefficients depend on time. 
Replacing the above expression of \^(t)) in eq.(l), we obtain the following set 
of equations for the amplitudes a m (t) 



= -iy2(n\V{t)\m)a m (t)exp(-iuj nm t) , (3) 
dt V 

where u nm = e m — e n are the Bohr frequencies. Combining eqs.(2) and (3), we 
find, 

da n (t) = Q (4) 



if n N and n ^ j; and 

da n {t) i 
dt ~ v^v' 



dt 



[1 - 5 nj ) aj(t) exp [+i (u jn + u sj ) t] 



+5 nj J2 a m(t) exp [-i (u jm + u sj ) t\) , (5) 

meN 
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if n G N or n = j. 



Before solving numerically the set of eqs.(4) and (5), it is important to un- 
derstand the qualitative behavior of these equations. Note that there are two 
time scales involved, a fast scale associated to the Bohr frequencies, and a 
slow scale associated to the amplitudes a n (t). Integrating the previous equa- 
tions in a time interval greater than the fast scale, the most important terms 
are those that have a very small phase, as the others average to zero. In this 
approximation the previous set of equations becomes 



dcij (t) 

dt ~ y/N 
da s (t) i 

dt ~ ~~ v^V 
da n (t) 



a s (t) , (6) 
a 3 (t), (7) 



dt 



~ for all n^s,j . (8) 



These equations represent two oscillators that are coupled so that their popu- 
lation probabilities alternate in time. As we notice the coupling is established 
between the initial and the searched for state. Solving those equations with 
initial conditions Oj(0) = 1, a s (0) = we obtain 

Pj ~ cos 2 (fi t) , (9) 
P s ~ sin 2 (fi t) , (10) 

where VL = -j=. It is important to notice that this approach is valid only if 
all the Bohr frequencies verify uj nm ^> Q. If we let the system evolve during 
a time r = | y/N , and we make a measurement immediately after that, the 
probability to obtain the searched state is one. For the case that the previous 
approximations are true, our method behaves qualitatively like Grover's, but 
the time r is twice the search time of Grover's algorithm [2]. This difference 
is not relevant since the resonance potential eq.(2) can be renormalized with 
a constant factor Vq as long as it does not invalidate the previous approxima- 
tions, uj nm ^> fi, where now fl — Vq/ \/N . 



3 Numerical results 



We have integrated numerically eqs.(4) and (5) for two Hamiltonians H , 
namely for the quantum harmonic oscillator and for the quantum rotor in 2D. 
The eigenvalues for these Hamiltonians are e m = e (m + |) and e m = s^rv? 
respectively. In both cases we take Eq — 1 and initial conditions aj(0) = 1, 
a n(0) = 0, for all n ^ j. The calculations were performed using a standard 
fourth order Runge-Kutta algorithm. Choosing an arbitrary eigenvalue for the 
energy of the searched state, we follow the dynamics of the set N. We have 
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verified, for several values of N , that the most important coupling is between 
the initial and the searched state, other couplings being negligible, as discussed 
in the previous section. 
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Fig. 1. Time evolution of the probability distribution of the initial and the search 
set states, for the quantum rotor, for N = {2, 10}. The initial state was taken to 
be j = 1, and the searched state s = 8. The full line shows the probability of the 
initial and the searched state. For all the others states of the set the probabilities 
take negligible values. 



The probability distribution of the search set and the initially loaded state are 
shown for the rotor case in Fig.l. Each panel shows the probability distribution 
at different times. It is clear that the total probability flows between the 
initial and the searched states. There is a characteristic time at which the 
probability of the searched state is maximum and its value is very near one. 
This time agrees with our theoretical prediction, r. Fig. 2 shows the oscillation 
of the probability flux between the initial and the searched states as a function 
of time for the quantum rotor. The temporal evolution was normalized for 
the characteristic time r. Fig. 3 shows a similar calculation for the harmonic 
oscillator. Here we have considered three sizes for the data set. Note two 
remarkable differences between Figs. 2 and 3: for the rotor the maximum values 
of the probabilities are one, while for the harmonic oscillator they are less than 
one, and decrease with time. Furthermore, while for the rotor the maxima 
are located in agreement with the theoretical prediction, for the harmonic 
oscillator the maxima show an increasing shift to the right of the predicted 
values as the search algorithm evolves. 

The differences between the two Hamiltonians mentioned above are due to the 
properties of their respective energy spectra: while the energy eigenvalues of 
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the rotor increase quadratically with the quantum number, for the harmonic 
oscillator they increase linearly. Consequently, the condition Q << uj nm , is 
better satisfied for the spectrum of the rotor than for the harmonic oscillator. 
As a consequence the agreement with the theoretical approximation is worse in 
this last case. However, as VI = ^= , if the size of the data set is increased, the 
agreement for the harmonic oscillator with the behavior predicted by eq.(10) 
is much improved, as shown in Fig. 3. Increasing N leads to a decrease in the 
relative distance between the energy spectrum values of Hq, allowing for a bet- 
ter resolution. In general, for a given Hamiltonian Hq with a discrete energy 
spectrum, the improvement obtained increasing N becomes negligible after 
some saturation value. Above this saturation value the algorithm presented 
here, has the same characteristics as Grover's algorithm. 




Fig. 2. Probability distribution for the initial (dashed line) and the searched (full 
line) states as a function of time for the quantum rotor Hq in the case of a data set 
of size N = 20 



The proportionality between the characteristic time r and y/~N for the rotor 
and the harmonic oscillator, is verified in Fig. 4. 

Taking iV greater than its saturation value, one can test the robustness of 
this search algorithm with respect to imprecisions in the eigenvalue of the 
searched state. Let us replace u S j by u S j = u S j + 5 in eqs.(5), where we take 
5 smaller than Q . 
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Fig. 3. Probability distribution for the harmonic oscillator Hq as a function of time. 
The dashed line corresponds to the initial state and the full line to the searched 
state. The size of the searched set are: (a) N = 20, (b) N = 60 , (c) N = 100. 
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Fig. 4. Time r at which the probability of the searched state is maximum as a 
function of the square root of the dimension of the search set. Squares correspond 
to the harmonic oscillator, and circles to the rotor. 

Within the this approximation the probability of the searched state evaluated 
in t — t is, 



P S (N,5) = 




We define the resonance width A, as the value of 5 for which the probability 
of the searched state falls to half of the value for 5 = 0. 
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Fig. 5. Probability to find the searched state as a function of the energy error 5 , 
for the rotor , with N = 10 (dashed line) and N = 100 (full line). A is the width of 
the resonance. 



The calculation of the probability of the searched state as a function of 5 is 
presented in Fig. 5. It can be observed that the curves are in agreement with 
eq.(ll), showing a symmetrical behavior about S = and a sharpness that 
depends on N. For small 5 the probability of the searched state remains large 
enough. Then, we deduce that a small error in u S j does not affect drastically 
the search algorithm. 



1/Vn 



Fig. 6. Resonance width as a function oil/y/N. The points are results of the numer- 
ical calculation. The dashed line corresponds to the linear regression of the data. 



The dependence of the resonance width with N is presented in Fig. 6. The 
numerical calculation fits with a hyperbola defined by Ay/N ~ 1.599. This 
result is in agreement with eq.(ll) when P s (JV, 8) 



~ \ with 5 = A. 
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4 Connection with a discrete time search algorithm 



We now build a discrete version of the continuous time algorithm developed in 
the previous sections. Note that the time r can be expressed by r = n/ (2ft) = 
VNuj S jT/4, where T = 27t/uj S j is the period of interaction potential, eq.(2). 
Then one can use Floquet's theory [10] to obtain the matrix of the Floquet 
evolution operator Uf — Uf(T). If one defines an unitary evolution operator 

Ud = (Uf) 4 ) the searched state is obtained by applying XJ D O(vN) times to 
the initial state. If uj S j/A is not an integer number, we round it to the nearest 
integer in the definition of Ud- 

To construct the Floquet operator, we must look for the amplitudes of the 
Floquet states that are also solutions of the eqs.(5). These amplitudes have 
the form 

a m (t) = exp (-i(\ - e m )t) (m\ <j> x (t)) , (12) 

where \(j>\(t)) = \<j>\(t + T)) is the X th Floquet state and A is the corresponding 
Floquet parameter. As \<j>\(t)) is periodic, it can be expanded as a Fourier series 

oo 

(m\(f>\(t))= A k (m, \) exp (ilusjt) . (13) 

k=— oo 

Using the eqs (12,13) in eq.(5) we obtain set of equations for the parameter A 
and for the amplitudes A k 



1 N 

(-A + ej + kuj sj ) A k (j, A) = — = £ A k+1 (m, A), (14) 

ViV m=1 

(-A + e n + ku sj ) A k (n, A) = y^A^j, A), for all n E N. (15) 

v N 

Replacing A k+ i from eq. (15) into eq. (14), the characteristic equation for A 
is obtained 



N 1 

xN= ^ (16) 

m=l x ' ^sm 

where x = —A + Ej + kuj S j. This equation has + 1 solutions for both x and 
A. Note that the indetermination in A introduced by the term kuj S j is superflu- 
ous since the Floquet evolution operator is evaluated at t — T and therefore 
exp(— ikujgjT) = 1. Once we have the values for A we obtain the Fourier coef- 
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ficients Ak(m, X) and, using eq.(13), (m\ 4>\(t)). 



The Floquet matrix has the form 



(Uf) 



mn 



(T) 



^exp(-iAT) (n|0 A (O)> (0 A (O)|m>, 



(17) 



A 



and, therefore, 



(^) 



mn 



(To) 



^exp(-ar )(n| A (O))(0 A (O)|m>, 



(18) 



A 



where To = u S jT/4. Interpreting Ud as a rotation operation with angle XTq, 
its application to the initial state 0(y/~N) times, maximizes the probability 
to find the searched state. Then for any H one can build a discrete time 
algorithm to perform the search. 



5 Conclusions 

We have developed a new insight in generating a continuous time quantum 
search algorithm using a characteristic of quantum mechanics, quantum res- 
onances. This algorithm behaves like Grover's algorithm; in particular the 
optimal search time is proportional to the square root of the size of the search 
set, y/~N, and the probability to find the searched state oscillates periodically 
in time. The efficiency of this algorithm depends on the spectral density of 
the Hamiltonian H . A bigger separation between the energy levels maximizes 
the probability of the searched state and allows for a better precision. For 
any Hamiltonian with discrete energy spectrum, independently of its spectral 
density, this algorithm can be implemented taking a large enough search set. 

The algorithm was shown to be robust when the energy of the searched state 
has some imprecision. However the improvement in the accuracy of the search 
by increasing N, and the error margin are bounded by the relation Ay/N = 
constant. This means that for a large N it demands a good precision in the 
eigenvalue of the searched state. Noting that A is an energy variation, and 
y/N is the time needed for the measurement, the previous relation is simply 
the Heisenberg uncertainty principle. 

We have found a simple way to build a discrete time algorithm on the basis of 
our continuous time search algorithm, which strongly suggests the equivalence 
between search algorithms with discrete and continuous time. 

It has been recently shown that Grover's algorithm may not be directly ap- 
plicable to search in a physical database [7]. Indeed, it has also been pointed 
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out that it would be neither technologically nor economically reasonable to 
build database search engines basen on Grover's search algorithm [11]. Similar 
questions might be raised about the algorithms presented here. 

Regarding the implementation of these algorithms, we note that the quan- 
tum kicked rotor has been experimentally realized using ultra-cold atom traps 
and some experiments have focused on the resonant case [12]. Furthermore, 
we have recently shown that the discrete quantum random walk on the line 
has the same dynamics as the kicked rotor in resonance [13,14]. Several sys- 
tems have been proposed as candidates to implement quantum random walks 
[15,16,17,18,19]. When these devices are constructed, they may be employed 
to simulate the search algorithms proposed here. 
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support from the Brazilian National Research Council (CNPq) and FAPERJ 
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